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This  paper  presents  an  analysis  of  the  fuel  flexibility  of  a  methane-based  solid  oxide  fuel  cell-gas  tur¬ 
bine  (SOFC-GT)  hybrid  system.  The  simulation  models  of  the  system  are  mathematically  defined.  Special 
attention  is  paid  to  the  development  of  an  SOFC  thermodynamic  model  that  allows  for  the  calculation 
of  radial  temperature  gradients.  Based  on  the  simulation  model,  the  new  design  point  of  system  for  new 
fuels  is  defined  first;  the  steady-state  performance  of  the  system  fed  by  different  fuels  is  then  discussed. 
When  the  hybrid  system  operates  with  hydrogen,  the  net  power  output  at  the  new  design  point  will 
decrease  to  70%  of  the  methane,  while  the  design  net  efficiency  will  decrease  to  55%.  Similar  to  hydrogen, 
the  net  output  power  of  the  ethanol-fueled  system  will  decrease  to  88%  of  the  methane  value  due  to 
the  lower  cooling  effect  of  steam  reforming.  However,  the  net  efficiency  can  remain  at  61%  at  high  level 
due  to  increased  heat  recuperation  from  exhaust  gas.  To  increase  the  power  output  of  the  hybrid  system 
operating  with  non-design  fuels  without  changing  the  system  configuration,  three  different  measures  are 
introduced  and  investigated  in  this  paper.  The  introduced  measures  can  increase  the  system  net  power 
output  operating  with  hydrogen  to  94%  of  the  original  value  at  the  cost  of  a  lower  efficiency  of  45%. 
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1.  Introduction 

High  temperature  solid  oxide  fuel  cells  (SOFCs)  seem  to  be  the 
most  efficient  device  for  direct  conversion  of  fuel  chemical  energy 
into  electricity  [1].  Its  efficiency  can  further  be  increased  by  inte¬ 
grating  it  with  a  gas  turbine  (GT).  This  kind  of  hybrid  SOFC-GT 
power  plant  is  an  attractive  near-term  option,  as  it  can  allow  effi¬ 
ciencies  of  over  60%,  even  for  small  power  outputs  (200-400  kW) 
[2]. 

In  most  published  papers,  methane  is  the  primary  fuel  for  SOFC- 
GT  hybrid  systems  [3-6].  Methane  is  the  main  component  of  natural 
gas,  which  is  widely  available  in  densely  populated  areas.  However, 
it  is  a  fossil  fuel  with  limited  resources,  and  cannot  be  produced 
from  renewable  sources.  Thus,  the  importance  of  renewable  or 
synthetic  fuels  is  expected  to  improve  during  coming  decades. 

Another  issue  with  natural  gas  is  that  it  might  not  be  economical 
to  deliver  gaseous  fuels  to  sparsely  populated,  remote  areas.  Liq¬ 
uid  fuels  are  advantageous  due  to  their  higher  volumetric  energy 
density,  which  facilitates  trailer  distribution  and  on-site  storage. 
Due  to  their  scalability  and  efficiency,  fuel-cell  power  plants  are 
considered  an  interesting  option  for  power  supply  in  sparsely  pop¬ 
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ulated,  remote  areas.  In  this  paper,  two  alternative  fuels  are  chosen 
to  fuel  the  hybrid  system,  namely  hydrogen  and  ethanol.  Hydrogen 
is  widely  expected  to  become  an  important  energy  carrier  in  the 
future  because  it  can  be  produced  from  almost  any  primary  energy 
source  with  acceptable  efficiency.  Ethanol  is  a  liquid  fuel  with  high 
volumetric  energy  density  and  can  be  produced  from  biomass. 

Thus,  it  is  important  to  determine  whether  an  SOFC-GT  hybrid 
system  designed  for  methane  fuel  can  be  easily  modified  to  other 
fuel  sources  once  the  focus  turns  from  natural  gas  to  other  renew¬ 
able,  novel  gases.  In  recent  years,  many  researchers  have  studied 
the  feasibility  of  SOFC  operating  on  biogas  or  other  suitable  fuels 
such  as  ethanol  and  methanol  [7-9].  Most  hold  the  opinion  that 
such  fuels  can  be  directly  fed  and  internally  reformed  in  an  SOFC 
stack.  Laboratory  tests  performed  by  Jenner  et  al.  confirmed  that  the 
efficiency  of  SOFC  drops  by  approximately  5%  when  the  mole  frac¬ 
tion  of  methane  in  the  fuel  decreases  from  70  to  30%  [10].  Yi  et  al. 
showed  that  the  electrical  efficiency  of  an  integrated  SOFC  drops 
around  1.1%  if  biogas  is  used  instead  of  natural  gas  [11].  Bo  Huang 
et  al.  fabricated  a  single  cell  of  SOFC  and  tested  its  electrochemical 
performance  fed  by  ethanol  steam  [12].  Few  papers  have  focused 
on  the  necessary  modification  and  performance  of  SOFC-GT  hybrid 
system  operating  with  non-designed  fuels. 

This  paper  discusses  the  performance  of  a  methane-based  SOFC- 
GT  hybrid  system  operating  with  non-designed  fuels.  Detailed 
simulation  models  of  all  components  of  the  system  are  mathemat- 
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Nomenclature 

Ae i  active  surface  of  fuel  cell  (m2) 

A  surface 

CpR  specific  internal  energy  (kj  kg-1 ) 

c  heat  capacity  (Jkg-1 1<)  or  mole  concentrations  of 

species 

F  Faraday  constant  (C)  or  molar  flow  rate  (mol  s-1 ) 

G  gas  flow  rate  (kgs-1) 

h  heat  transfer  coefficient  (Wm-2I<-1)  or  specific 

enthalpy  (kj  kg-1) 

Ah  enthalpy  change  of  the  reaction  (J  mol-1 ) 

I<  equilibrium  constant  of  reaction 

k  thermal  conductivity(  W  K-1  m-1 ) 

I  current  density  (A  m-2 ) 

l  tube  length 

fn  flow  rate 

P  pressure  (Pa)  or  power  density  ( W  m-2 ) 

p  partial  pressure  (Pa) 

q  radiation  heat  flow  (W) 

R  gas  constant  (J  mol-1  K-1 )  or  resistance  (£2) 

r  radius  (m) 

r  reaction  rate  (mol  s-1 ) 

T  temperature  (I<) 

U  voltage  (V) 

u  flow  velocities  (ms-1) 

W  power 

x  mole  fraction 

Greek  symbols 
q  efficiency 

p  density  (kg  m-3 )  or  ohmic  resistance  ( £2  m) 

e  compression  pressure  ratio 

8  expansion  pressure  ratio  or  thickness  (m) 

Subscripts 

ADT  air  delivery  tube 

c  compressor 

Fu  fuel 

/  fuel  channel 

MEA  membrane  electrolyte  assembly 

op  over  potential 

ocp  open  circuit  potential 

Ref  reforming  reaction 

Shi  water-gas  shift  reaction 

T  turbine 


ically  defined,  paying  more  attention  to  the  SOFC  thermodynamic 
model  developed.  Based  on  the  simulation  models,  the  feasibility 
and  steady-state  performance  of  a  methane-based  SOFC-GT  hybrid 
system  operating  with  hydrogen  and  ethanol  are  discussed.  Three 
different  methods  are  introduced  to  increase  the  power  output  of 
a  hybrid  system  operating  with  non-designed  fuel. 

2.  Mathematical  models 

2.1.  System  configurations 

The  SOFC-GT  hybrid  system  configuration  is  shown  in  Fig.  1. 
The  steam  required  to  support  reforming  reactions  (reactions  (I) 
and  (II)  in  Table  1 )  is  derived  from  the  anode  outlet  stream.  During 
the  start-up  procedure,  the  plant  requires  an  external  boiler  that 
produces  steam  for  the  reforming  reaction.  In  this  paper,  the  startup 
procedure  is  not  taken  into  account;  thus,  the  external  boiler  is  not 


Auxfuel  valve 


Fig.  1.  Configurations  of  SOFC-GT  hybrid  system. 


considered.  The  hybrid  system  mainly  includes  an  SOFC  subsystem, 
a  gas  turbine  subsystem,  two  exchangers,  a  combustor,  and  a  pre¬ 
reformer.  The  working  principle  of  the  plant  can  be  summarized  as 
follows. 

The  compressed  air  is  preheated  through  two  heat  exchangers 
(HE1  and  FIE2)  before  entering  the  cathode  of  the  SOFC  stack.  In 
the  anode  compartment,  fuel  is  mixed  with  anode  re-circulated 
stream  before  entering  the  pre-former.  After  being  preheated  and 
partly  reformed  in  the  pre-reformer,  it  is  carried  to  the  external  sur¬ 
face  of  the  SOFC  tube,  where  internal  reformation  occurs,  releasing 
hydrogen,  which  is  brought  to  three-phase  boundaries  (TPBs)  [13]. 
Preheated  air  is  fed  to  the  top  of  the  stack,  where  is  further  pre¬ 
heated  and  finally  carried,  by  air  delivery  tube  (ATD),  into  the  SOFC 
tube  as  shown  in  Fig.  2.  There,  it  participates  in  the  cathode  electro¬ 
chemical  semi-reaction.  The  electrochemical  reactions  occur  at  the 
TPB  of  both  electrodes,  producing  an  ionic  flow  through  the  elec¬ 
trolyte  and  electron  across  the  electrodes.  Electrical  energy  is  hence 
produced  associated  with  heat  generation  during  the  process.  The 
heat  generated  is  partly  used  to  reform  the  fuels,  partly  dissipated 
to  the  environment,  and  partly  used  to  heat  up  the  feedstock  gases. 
The  excess  air  and  unreacted  fuels  are  burnt  within  the  combustor 


Table  1 

Reactions  in  SOFC  tube. 


Steam  reforming  reaction 

CH4  +  H20  3H2  +  CO 

(I) 

C2H5OH  +  3H20  2C02  +  6H2 

(II) 

Water  gas  shift  reaction 

co+h2o^co2+h2 

(HI) 

Oxidation 

H2+02-  ->  H20+2e- 

(IV) 

Reduction 

02+4e-^202- 

(V) 
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Fig.  2.  Geometric  profile  of  SOFC  tube. 


to  increase  stream  temperature  partly.  HE2  is  placed  before  the  GT 
inlet  to  preheat  the  SOFC  inlet  air.  The  high  temperature-pressured 
gas  then  expands  in  the  turbine  and  produces  mechanical  work 
to  drive  the  air  compressors  and  the  directly  coupled  electrical 
generator. 

2.2.  SOFC  model 

The  SOFC  model  developed  in  this  study  is  based  on  a  tubu¬ 
lar  design.  The  geometry  of  the  SOFC  is  based  on  the  design  of 
Siemens-Westinghouse  [14,15];  the  geometric  profile  is  shown  in 
Fig.  2.  The  model  developed  in  this  paper  is  based  on  the  following 
assumptions. 

1.  The  model  described  a  bundle  of  non-interacting  fuel  cell  tubes, 
no  stack  effects  are  considered. 

2.  The  gas  flows  are  modeled  as  plug  flows. 

3.  Laminar  flow  is  assumed  in  heat  transfer  models. 

4.  The  shift  reaction  is  assumed  to  be  always  at  equilibrium. 

5.  All  heat  sources  and  sinks  derived  from  the  reactions  attack  at 
the  anode  surface. 

The  SOFC  model  presented  in  this  paper  consists  of  both  electro¬ 
chemical  and  thermal  models.  The  thermodynamic  model  includes 
convectional  heat  transfer  between  gases  and  solids  as  well  as 
radiation  and  conduction  between  solid  materials  is  developed, 
allowing  for  the  calculation  of  radial  temperature  gradients.  Kinet¬ 
ics  of  reforming  and  electrochemical  reactions  is  included  in  the 
electrochemical  model. 

2.2.1.  Electrochemical  reactions 

In  an  SOFC,  the  electric  power  is  generated  through  the  following 
reactions  as  shown  in  Table  1. 

The  reaction  rate  of  the  CH4  reforming  reaction  (I)  is  determined 
by: 

rRef( I)(Z)  =  fcact  '  e-(£ac.)/(K-WrE,m,z))  .  Xfu(z,  CH4)  •  PFu,ref  •  A*  (1 ) 

where  Fact  is  the  activation  energy  while  /<act  is  the  pre-exponential 
constant. 

The  value  of  Fact  and  /<act  have  been  given  in  Ref.  [16] 
Fact  =  82  x  103  J  mol-1;  kact  =  0.04274  mols”1  m~2  Pa-1. 


Such  formation  is  considered  typically  for  the  DIR-SOFC 
performance. 

C2H5OH  reforming  reaction  is  much  more  complex  than  that 
of  CH4  which  has  been  described  in  many  papers  such  as  [17]. 
The  ethanol  reforming  reaction  kinetic  model  derived  by  Abayomic 
Akande  are  adopted  in  this  paper.  The  model  is  of  the  form 

fjwrofcO  =  feo  ■  Xn(z,  C2H5OH)  •  e— (Eo)/(R- WrE,m,z))  pFu  ref .  Ael  (2) 

where  k0  is  the  pre-exponential  constant;  F0  is  the  activation 
energy;  n  is  the  order  of  reaction  with  respect  to  ethanol.  The  value 
of  Fact,  n  and  kact  have  been  given  based  on  experimental  data  in 
Ref.  [18],  in  which  k0  =  3A2x  10~2;  n  =  0.43;  F0  =  4.41  x  103. 

The  reaction  rate  of  the  shift  reaction  can  be  calculated  by: 

rshi(z)  .  104  (*Fu<z,  CO)  xF„(z.  CO)  -  C°4 1 )  : 

X  z  e  [0, 1)  (3) 

where  Kshi(z)  is  the  equilibrium  constant  and  can  be  calculated  by: 
(4)ln(Kshi(z))  = 

The  electrochemical  reaction  rate  is  calculated  from  the  dis¬ 
tributed  current  density. 

nacU)  =  (5) 

2.2.2.  Heat  transfer  and  radiation 

Both  fuel  reforming  and  water-gas  shift  reactions  are  endother¬ 
mic,  and  the  required  heat  of  the  reactions  is  calculated  by: 

A/iRefdjCz)  =  3  •  (o  +  h{ H2))  +  (h°(C 0)  +  Ji(CO)) 

-  (h°(CH4)  +  h(CH4))  -  (Ji°(H20)  +  Ji(H20))  (6) 

A/lReKH P)  =  4  .  (0  +  /i(H2))  +  2  ■  (h°( CO)  +  h(CO))  -  (h°(C2H5OH) 
+  h(C2H50H))-(/l°(H20)  +  /l(H20))  (7) 

Abshi(z)  =  3  ■  (0  +  h( H2))  +  (h»( CO)  +  Ji(CO))  -  (h°(CO)  +  h(CO)) 
-(fl°(H20)  +  h(H20))  (8) 
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For  the  electrochemical  reaction,  part  of  the  enthalpy  change  is 
turned  into  electric  power. 

A hBc(z)  =  (h°(H20)  +  h(H20))  -  i(h?( 02)  +  h(02))  -  (0  +  h( H2)) 
+  2  UF  (9) 


The  boundary  heat  transfer  between  the  ADT  and  ADT  air  is 
calculated  by: 

^ADT  •  — AD9rApjT,i  ’  ^  =  ^c,ADTAir-ADl(z)  *  (TaDt(^ADT,i?  z) 

-TADTAirW);  Z  G  (0,  /)  (10) 


The  heat  conduction  in  the  radial  and  axial  direction  inside  the 
ADT  is: 

„  dIADT(r,  z)  ( d2TAm{r,z)  3(r  •  dTAm(r,z))\  # 

CADT  •  PADT - =  fcADT  •  ^  9z2  +  r  .  ^  )  ' 

x  r  G  (rADT,i>  rADT,o)>  Z  G  (0, /)  (11) 


The  boundary  heat  transfer  between  the  anode  and  the  fuel 
(including  heat  of  the  reaction  and  radiation  to  the  reformer  and 
the  casing)  is  calculated  by: 

h  _  3TMEA(r0,  z) 

^MEA.r  '  K~ 

WMEA 

=  -hc,Fu-A(z)  •  (^MEA (f0,  Z)  -  Tpu(z)) 

^Ref^^Ref  +  ^Shi^^Shi  +  ^Elc^^Elc  +  9 rad, Ref  (z)  +  9rad,Cas(z) 
27tr0l 

(17) 

where  qrad,Ref(z)  and  <7rad,cas(z) are  the  radiation  heat  flows  from  the 
anode  surface  to  the  reformer  and  the  anode  surface  to  the  stack 
casing,  respectively: 

<?Ref(z)  =  TRef  •  '  (7e,Ref(z)4  “  T'r.Rerfz)4)  (18) 

^Rad.casC^)  =  ^Cas  ‘  ‘  (te,Cas(z)  —  l'r.Cas(^)  )  (19) 

The  heat  transport  in  the  fuel  is  calculated  by: 


5>u(z,  i)  ■  Cp,Fu(z)  •  r  dr^'Z)  +  UfuCD  ■ 
1=1  ' 


3TFu(z) 

9zf 


The  boundary  heat  transfer  between  ADT  and  cathode  air 
(including  radiation  from  the  cathode  surface)  is  calculated  by: 

^ADT  •  — AD^  ADT’n  ^  =  -hc,Air-ADl(z)  *  (7aDt(*ADT,o>  z)  -  TAir(z)) 
°,rADT 

+  qczmIiz)  (12) 

2tt  rADT)0/ 


The  heat  transport  in  the  cathode  air  is  calculated  by: 


^cAir(z,  i)  ■  CpiAir(z)  •  ( 
1=1  ' 


2  •  rADT)0 


dTAir  (z)  3TAir(z) 

fa  +  UAir(0  jgj— 


'  l,c,Air-ADl(z)  •  (lADl(rADT,o>  Z)  -  TAir(z)) 


+  “2 — -y - hc,Air-c(z)  •  (TMEA(rijZ)  -  TAir(z)), 

r,  -  rADT,o 

X  ze(0,  /)  (13) 

The  boundary  heat  transfer  between  cathode  air  and  cathode 
(including  ratio  to  the  ADT  surface)  is  calculated  by: 

^MEA.r  •  — ^A^  1  ’  ^  =  -hc,Air-c(z)  *  (^MEA (ri>  z)  “  ^Air(z)) 

wmea 

,  <7c-adt(z)  /  'i  a  \ 


=  ^  •  hc,Fu-A(z)  •  (TMEA(r0,  Z)  -  Tpu(z)), 

^Fu 

X  zg  (0,  /)  (20) 

Assuming  r0  -  r,-  <<c  rit  the  mean  temperature  over  the  MEA  can 
be  calculated  by. 


1-lTo-n ) 


TMEA(r,  z)drdz 


The  local  heat  gradient  maxima  in  the  radial  of  the  MEA  were 
calculated  by: 


rMEA,max  —  Hiax 


3TMea(^  z) 


2.2.3.  Mass  transfer  model 

On  the  air  side,  there  is  only  one  sink  term  exit  for  oxygen,  and 
no  source  or  sink  terms  exit  for  other  components. 

dcAir(z,02)  ,  „  dcAir(z,  02) 

— dt —  +  “Air(,) — aZ — 


(~l/2)rEic(z) 

/  ■  (r?  —  r?nT  )n 


V  z  e  [0,  /) 


dcAir(z,  i)  +  UAir(/)9cAg,  I)  =  0  ze[0,/),  i  e  {N2,  H2,  CH4,  CO, 

H20,C02,C2H50H)  (24 


where  <7c-adt(z)  is  the  radiation  heat  flow  between  the  cathode  and 
the  ADT,  which  can  be  calculated  by: 

<7c-adt(z)  =  ^c-adt  •  •  (Tmea (n,  z)4  -  TAm{rAm>0,  z)4)  (15) 

The  heat  transfer  in  the  radial  and  axial  directions  inside  the 
MEA  (note  different  heat  conductivities  in  the  radial  and  axial  direc¬ 
tions)  is  calculated  by: 

drMEA(r,  z)  32rMEA(r,  z)  3  (r  •  3rADT(r,  z)) 

dt  =  W  ■  — -  +  fcMEA,r  •  r  . 

X  r  g  (r,-,  r0),  zg(0, /)  (16) 


On  the  fuel  side,  there  is  no  source  or  sink  terms  exit  for  nitrogen 
or  oxygen. 

dCFaJf°+uFa(ndCF^’,)=0;  ze[0,  /),  ie{ N2;02)  (25) 

The  mass  transfer  equations  for  other  species  are: 

dcEu(z,  H2)  3cFu(z,H2) 

- dt -  +  UFu(,)  3zc 

_  2rRef(i)(z)  +  4rRef(I)(z)  +  rshi(z)  -  rElc(z)  . 
l  ■  AEu 

x  z  g  [0,  /)  (26) 
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dcFu(z,  CH4)  9Cpu(2,  CH4)  ^Ref(I)(2).  rn  n 

- dt - +UFu(0 - dzc - =  lAfu  •  Ze[0J) 

dCFu(Z,  C2H5OH)  dcFu{Z,  C2H5OH4)  -^Ref(II)(z) 

- dt -  +  UFu(,) - fee - =  l.Ato  : 

x  ze[0, 1) 


(27) 


(28) 


dcFu(z,  H20) 
dt 


+ 

OZ  c 


— rRef(i)(z)  -  rRef(n)(z)  -  rshi(z)  +  rEic(z) . 
I  ■  AFu 


x  ze[0,  /) 


(29) 


dcFu(z,  CO) 
dt 


+  UFu(0 


9cFu(z,  CO) 
9zc 


^Ref(I)(z)  +  2rRef(ii)(z)  -  fshi(z)  . 
I  •  AFu 


XZG  [0,  /) 


(30) 


dcFu(z,  C02) 
dt 


+  wFu(/) 


9cFu(z,  C02) 
dzc 


fshi(z). 
I  •  AFu 


ze[0J) 


(31) 


2.2.4.  Electrochemical  model 

For  SOFC,  electrochemical  models  are  coupled  with  heat  and 
mass  transfer  models  to  determine  the  distribution  of  local  current 
density  and  electrode  over  potentials  inside  the  SOFC. 

The  relationship  between  cell  current  density  and  voltage  is 
written  as  follows: 

j  =  _ tVp  ~  Id  op _ 

^ohm  +  ^anode  +  ^cathode 

where  Li0cp  is  the  function  of  SOFC  operation  temperature  and  gas 
partial  pressure  in  anode  and  cathode,  which  can  be  regulated  as 
follows: 


U  U°  I  RT  Inf  Pf’H2(Pa'°2) 

l/ocp-l/H2  +  -lnl  — 


(33) 


where 

U^2  =  1.2723  -2.7645e-4  x  T  (34) 

U^2  is  the  ideal  voltage  for  hydrogen  oxidization  [19],  pstd  is  the 
standard  pressure. 

The  output  power  density  of  SOFC  can  be  calculated  by: 


Psofc  =  V  x  /  (35) 

As  both  electrodes  are  normally  good  conductors,  a  constant  cell 
operation  voltage  Uop  throughout  the  cell  is  normally  considered 
[20]. 

Ohmic  losses  are  caused  by  resistance  to  conduction  of  ions  and 
electrons,  and  contact  resistance  between  cell  components. 


^ohm  —  Pe^e  +  Pc$  c  +  Pa^a  (36) 

The  activation  polarization  terms  are  controlled  by  the  electrode 
reaction  kinetics  of  the  respective  electrodes.  They  represent  volt¬ 
age  loss  incurred  due  to  the  activation  necessary  for  charge  transfer. 
Over  potentials  at  the  cathode  and  anode  are  assumed  indepen¬ 
dent  of  the  local  current  density.  The  model  derived  by  Achenbach 
is  adopted  in  this  paper  [21  ]. 


Fig.  3.  CH4  fueled  SOFC  model  validation. 


Current  density  (Am'2) 

Fig.  4.  C2H5OH  fueled  SOFC  model  validation. 


2.2.5.  Model  validation 

The  SOFC  model  developed  in  this  paper  is  validated  using  the 
experimental  data  published  in  Ref.  [22].  The  input  parameters  and 
fuel  utilization  values  of  the  experimental  data  were  simulated  with 
the  model,  giving  voltage  and  power  as  results.  Input  parameters 
and  results  are  shown  in  Table  2.  The  simulation  results  match  both 
experimental  datal  and  experimental  data  2  with  errors  of  only 
1.01  and  1.09%,  respectively. 

In  order  to  further  validate  the  SOFC  sub-model,  the  polariza¬ 
tion  curves  generated  by  the  simulation  model  are  compared  with 
experimental  ones  for  different  fuels  as  shown  in  Figs.  3-5.  Results 
shows  that  the  developed  model  can  achieve  errors  lower  than  3%. 
Unfortunately,  we  have  not  acquired  the  SOFC  stack  experimental 
data,  so  the  validation  of  the  stack  performance  is  not  yet  possi¬ 
ble.  CH4  fueled  SOFC  sub-model  was  validated  by  experimental 
data  published  in  [23].  C2FI5OFI  fueled  SOFC  sub-model  was  val¬ 
idated  with  experimental  results  with  the  operating  temperature 
at  1 1 23 1<  and  1 073  K  published  in  [24].  H2  fueled  SOFC  sub-model 
was  validated  with  experimental  data  from  Siemens  published  in 
[25]. 

2.3.  Gas  turbine  model 

The  micro  gas  turbine  selected  in  this  paper  is  C30  from  CAP¬ 
STONE,  which  is  a  single-shaft  micro  gas  turbine  equipped  with 
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Table  2 

Validation  of  the  SOFC  model. 


Experimental  data  1 

Simulation  results 

Experimental  data  2 

Simulation  results 

Pressure  (bar) 

1.05 

1.05 

3.5 

3.5 

Fuel  flow  rate  (mol/s) 

1.511  x  10“3 

1.511  x  10“3 

2.287  x  10“3 

2.287  x  10“3 

Air  flow  rate  (mol/s) 

1.055  x  10-2 

1.055xl0-2 

1.29xl0-2 

1.29  x  10-2 

Inlet  fuel  temperature 

823  K 

823  K 

860  K 

860  K 

Inlet  air  temperature 

1 1 04  K 

1 1 04  K 

1048  K 

1048  K 

Fuel  utilization 

Inlet  fuel  composition  (mole  fraction) 

69% 

69% 

69% 

69% 

h2 

25.8% 

25.8% 

22.6% 

22.6% 

h2o 

28.4% 

28.4% 

33.4% 

33.4% 

ch4 

11% 

11% 

13.1% 

13.1% 

CO 

5.7% 

5.7% 

5.7% 

5.7% 

co2 

22.8% 

22.8% 

24.1% 

24.1% 

n2 

Results 

6.3% 

6.3% 

1.1% 

1.1% 

Voltage  (V) 

0.69 

0.697 

0.639 

0.632 

Power  (W) 

104.8 

105.86 

157 

155.28 

Error 

- 

1.01% 

- 

1.09% 

Current  density  (  A  m'2) 

Fig.  5.  H2  fueled  SOFC  model  validation. 


Table  4 

Thermodynamic  equations  for  compressor  and  turbine  models. 


Compressor  outlet 
pressure 

Compressor  outlet 
temperature 
Compressor 
consumed  power 
Turbine  outlet 
temperature 
Turbine  output 
power 


(37)Pout  =  ePin 


(38  )Tc,out  =  rc,in  +  Tc,in(x{k-m  -  1  )/!)c 
(39)WC  =  4.18CpGkrCj„(£(k-1)/fc  —  1  )tr)c 


(40)Tt.ollt 


=  rt,i„ 


( r_p*YT~mT 


r  /  p  \  (Xt— 1)/Xt-i 

(41)WC  =  4.18cp,RCTTtjin  [l  -  )  Jrj t 


Generator  efficiency  in  this  paper  is  assumed  constant.  The 
power  output  from  the  gas  turbine  can  be  calculated  by: 

WGT  =  rjgen(rjtWt  -  Wc)  (42) 

where  r]gen  is  the  generator  efficiency,  r]t  is  the  turbine  mechanical 
efficiency,  and  Wt,  Wc  are  turbine  output  power  and  compressor 
consumed  power,  respectively. 


a  centrifugal  compressor  and  a  radical  turbine.  Major  parameters 
of  C30  gas  turbine  are  summarized  in  Table  3.  The  compressor  and 
turbine  models  provide  the  standard  thermodynamic  equations  for 
the  input-output  relationship  of  the  working  fluid  based  on  pres¬ 
sure  ratio  and  isentropic  efficiency,  as  shown  in  Table  4.  The  latter 
values  are  determined  by  a  performance  map,  which  describes  the 
compressor  and  turbine  by  introducing  pressure  and  efficiency  as 
a  function  of  mass  flow  rate  and  shaft  speed.  Performance  maps  of 
the  compressor  are  shown  in  Fig.  6. 


2.4.  Combustor  model 

Assume  that  all  residual  combustible  compounds  such  as 
C2H5OH,  CH4,  CO,  and  H2  formed  in  the  fuel  cell  stack  could  be 
fully  oxidized  in  the  combustor.  The  associated  chemical  reactions 
are  as  follows: 

C2H5OH  +  302  =  2C02  +  3H20  +  AHrxn,CH4  (43) 

CH4  +  202  =  C02  +  2H20  +  AHrxn>CH4  (44) 


Table  3 

Parameters  of  the  C30  gas  turbine  at  design  point. 


Parameters 

Value 

Pressure  ratio 

3.2 

Compressor  isentropic  efficiency 

80% 

Air  flow  rate 

0.31  kgs-1 

Fuel  flow  rate 

0.0024  kg  s_1 

Turbine  inlet  temperature 

1173K 

Turbine  isentropic  efficiency 

90% 

Turbine  exit  temperature 

548  K 

Turbine  exit  pressure 

103  KPa 

Recuperator  effectiveness 

90% 

Electric  power  generated 

30  kW 

Electric  efficiency 

26% 

CO  +  io2  =  C02  +  AHrxnco  (45) 

H2  +  lo2=H20  +  AHrxn,H2  (46) 


The  temperature  of  the  product  gas  (Tp)  can  be  calculated  based 
on  a  simple  energy  balance  on  the  combustor  control  volume. 

+  Ahrxn,c2H5OH  +  Ahrxn,cH4  +  Afarxn  co  +  Ahrxn,H2 


(47) 
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Fig.  6.  Performance  maps  of  compressor. 


2.5.  Heat  exchange  model 

The  heat  exchanger  model  presented  in  this  paper  is  simulated 
based  on  the  s-NTU  method  [26],  implementing  temperature- 
dependent  specific  heat. 

An  iterative  procedure  is  implemented  to  assess  the  out¬ 
let  temperatures,  based  on  the  calculation  of  average  hot 
and  cold  fluid-specific  heats.  These  parameters  depend  on  the 
unknown  HE  outlet  temperature.  Thus,  estimated  values  for 
the  temperatures  are  selected,  allowing  for  the  calculation  of 
average  specific  heats  and  consequently,  through  the  s-NTU 
method,  the  outlet  temperature  is  re-calculated.  The  procedure  is 
stopped  when  the  convergence  criterion  on  such  temperature  is 
satisfied. 

3.  Results  and  discussion 

3.1.  System  performance 

3.1.1.  Performance  at  design  point 

The  hybrid  system  operation  parameters  at  design  points  fed 
by  methane,  hydrogen,  and  ethanol  are  shown  in  Table  5.  Fig.  7 


Table  5 

Design  parameters  of  hybrid  system. 


Parameters 

Methane 

Hydrogen 

Ethanol 

ft  fuel  (mols-1) 

0.2996 

0.767 

0.1708 

ftairtkgS”1) 

0.31 

0.307 

0.34 

Pnet  (kW) 

150 

105 

132 

/(Am-2) 

2980 

1910 

2532 

TIT  (K) 

1173 

1069 

1132 

Pmea  (K) 

1273 

1273 

1273 

Ysc 

2 

- 

1.7 

7TC 

3.2 

3.04 

3.13 

4  net 

62.29% 

56.59% 

60.5% 

shows  the  process  parameters  at  design  points  of  the  hybrid 
system  operating  with  all  three  fuels.  To  evaluate  the  perfor¬ 
mance  of  the  system  fed  by  new  fuel,  new  design  points  must 
be  defined  for  each  fuel  at  first.  The  gas  turbine  is  only  slightly 
influenced  by  the  fuel  choice,  thus,  the  original  air  flow  rate 
from  the  methane-fueled  system  is  chosen.  SOFCs  are  the  most 
important  component  in  the  system,  contributing  most  of  the 
net  power.  To  make  SOFCs  operate  at  their  design  point,  the  fuel 
flow  for  each  new  fuel  is  designed  with  the  objective  of  reach¬ 
ing  the  same  SOFC  operating  temperature  as  with  methane,  which 
is  1273  K. 

When  the  hybrid  system  is  operating  with  hydrogen,  the  design 
calculation  returns  a  fuel  flow  rate  of 0.767  mol  s-1 , 2.56  times  more 
than  that  of  methane.  Assuming  that  fuel  utilization  is  constant,  this 
value  corresponds  to  an  SOFC  current  density  of  only  1910  Am-2. 
The  reason  for  the  changed  relationship  between  SOFC  design  cur¬ 
rent  density  and  operating  temperature  is  the  lack  of  cooling  effect 
from  steam  reforming.  Higher  air  flow  rate  should  be  required  for 
the  hydrogen-fed  hybrid  system  to  reach  the  original  design  current 
density  at  an  acceptable  temperature  level.  The  reduction  current 
implies  a  reduction  of  design  power  to  0.7  of  the  value  reached 
with  methane.  The  net  efficiency  is  reduced  to  56.59%.  The  lower 
efficiency  is  due  to  the  higher  cooling  air  demand  per  power,  and 
thus  higher  heat  loss.  The  same  effect  has  been  observed  by  other 
researchers  [27]. 

The  design  calculation  for  ethanol  fuel  returns  a  relative  fuel 
flow  rate  of  57%  of  the  methane  value,  that  is,  0.1 708  mol  s-1 .  With 
higher  molar  hydrogen  yield  of  ethanol,  this  value  corresponds  to 
an  SOFC  current  density  of  2549  Am-2.  The  changed  relationship 
between  SOFC  temperature  and  current  density  is  caused  by  lower 
specific  heat  per  mole  of  hydrogen  required  by  the  reforming  reac¬ 
tion.  The  reduced  current  density  implies  a  reduction  of  net  power 
to  88%  of  the  value  reached  with  methane  fuel,  at  a  net  efficiency 
of  60.5%.  These  results  are  achieved  based  on  the  assumption  that 
ethanol  vaporization  is  accomplished  by  exhaust  gas  heat.  If  an  elec¬ 
tric  heater  is  used  for  this  purpose,  the  net  power  will  decrease  to 
80.02%  of  the  methane  value  and  the  net  efficiency  will  decrease  to 
55.4%. 

3. 1.2.  Part-load  performance 

In  the  SOFC-GT  hybrid  system,  SOFC  is  the  most  impor¬ 
tant  component  that  contributes  most  of  the  net  power.  To 
make  SOFC  operate  at  its  design  point,  a  part-load  strategy 
with  a  constant  SOFC  mean  MEA  temperature  is  pursued  in  this 
paper. 

Fig.  8  shows  the  major  parameters  and  system  performance 
of  the  SOFC-GT  hybrid  system  operating  with  methane  at  part¬ 
load  conditions.  The  lowest  electrical  power  of  21%  of  the  designed 
value  (31.5  kW)  can  be  achieved  by  reducing  the  air  and  fuel  flow 
rates  at  constant  SOFC  operating  temperature.  At  partial  load  oper¬ 
ation,  less  fuel  is  fed  into  the  SOFC  stack  while  the  fuel  utilization 
is  assumed  constant.  Consequently,  the  unreacted  fuels  in  SOFC 
exhaust  decreases,  determining  the  reduction  of  turbine  inlet  tern- 
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Ethanol 

T  (K) 

T  (K) 

T(K) 

P  (bar) 

P  (bar) 

P  (bar) 

288 

288 

288 

1.01 

1.01 

1.01 

Fig.  7.  Process  parameters  of  hybrid  system  at  design  point. 


perature.  The  variation  of  the  gas  turbine  air  flow  rate  and  the 
reduction  of  turbine  inlet  temperature  (TIT)  cause  the  reduction 
of  both  GT  output  power  and  efficiency.  On  the  other  hand,  an 
increase  in  cell  voltage  is  obtained  because  of  the  current  den¬ 
sity  reduction  caused  by  the  reduction  of  the  fuel  flow  rate.  The 
high  voltages  achieved  by  the  cell  determine  the  lower  efficiency 
defects  of  the  SOFC  stack.  These  effects  are  in  contrast,  which  cause 
the  system  efficiency  to  be  in  a  steady-state  ranging  from  59.6%  at 
minimum  load  to  a  maximum  of  65.8%  at  intermediate  load  (59% 
power),  and  back  to  62.6%  at  the  design  point.  For  the  SOFC  tube,  the 
minimum  MEA  temperature  has  a  sharp  bend  at  68%  power,  where 
the  maximum  MEA  temperature  is  lowest  at  1294K.  The  highest 
local  gradient  in  radial  direction  occurs  at  the  design  point  with 
797  K  nrr 1 ,  caused  by  the  higher  cooling  effect  of  the  fuel  reforming 
reaction. 

Fig.  9  shows  the  major  parameters  and  system  performance 
of  the  SOFC-GT  hybrid  system  operating  with  hydrogen  at  part¬ 


load  conditions.  The  electrochemical  power  of  the  hybrid  system 
can  be  produced  in  the  range  14-70%  of  the  value  reached  with 
methane,  that  is,  21  -1 05  kW.  Compared  with  the  operation  param¬ 
eters  for  the  methane-fueled  system,  the  strongest  difference  is 
the  maximum  radial  temperature  gradient,  which  for  a  certain 
power  output  (105  kW,  70%  of  the  design  power  with  methane), 
is  approximately  three  times  higher  for  the  methane  system.  For 
both  options,  the  maximum  radial  gradient  occurs  at  the  cell  bot¬ 
tom,  where  air  and  fuel  meet  at  either  side  of  the  MEA  with 
different  temperatures.  Due  to  the  absence  of  an  endothermic- 
reforming  reaction,  the  temperature  of  the  entering  fuel  is  higher 
than  if  methane  is  used  and  the  difference  to  the  air  is  thus 
lower.  The  turbine  inlet  temperature  is  lower  with  hydrogen 
fuel  caused  by  lower  SOFC  exhaust  temperature,  both  at  equal 
power  and  air  flow  rate  are  compared.  This  causes  the  gas  tur¬ 
bine  to  be  operating  on  a  different  line  with  a  methane-based 
system. 
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Fig.  8.  Electrochemical  performance  and  operation  parameters  of  SOFC-GT  hybrid  system  fed  by  methane. 


Fig.  10  shows  the  major  parameters  and  system  performance 
of  an  SOFC-GT  hybrid  system  operating  with  ethanol  at  part¬ 
load  conditions.  Net  power  output  is  in  the  range  of  19-88%  of 
power  reached  by  methane-based  system  at  design  point.  At  88% 
power,  the  net  efficiency  is  60.5%,  lower  than  the  value  at  88% 
power  by  methane  system,  that  is,  63.7%.  Similar  to  the  hydrogen- 
fed  system,  the  ethanol  fuel  also  slightly  reduced  the  thermal 
load  to  the  fuel  cell,  as  the  maximum  local  MEA  temperature  is 
lower,  and  the  minimum  MEA  temperature  is  higher  than  that  for 
the  methane  case.  At  88%  power,  the  maximum  radial  tempera¬ 
ture  gradient  is  approximately  40%  lower  than  that  of  methane, 
which  is  a  consequence  of  the  lower  heat  consumption  of  the 
reforming  reaction.  The  turbine  inlet  temperature  is  higher  than 
that  of  H2,  but  lower  than  that  of  methane  at  equal  net  power 
output. 

Fig.  1 1  shows  the  gas  turbine  operational  line  in  the  compres¬ 
sor  map.  Although  the  air  flow  rate  decreases  with  the  reduction 
of  system  output  power,  the  reduction  of  TIT,  which  causes  a 
pressure  ratio  decrease,  ensures  that  the  compressor  surge  will 
not  occur  during  the  system  off-design  operation  range  for  all 
fuels. 


3.2.  Measures  to  increase  power  output  for  non- designed  fuel 

To  increase  the  power  output  of  the  hybrid  system  operating 
with  non-designed  fuel,  the  following  measures  were  introduced 
and  investigated  in  steady-state. 

•  Case  1 :  Using  the  auxiliary  fuel  to  improve  the  TIT  to  its  design 
value. 

•  Case  2:  Exhaust  gas  bypassing  the  HE1  to  decrease  heat  recuper¬ 
ation. 

•  Case  3:  Exhaust  gas  bypassing  the  HE1  and  using  auxiliary 
fuel. 

As  discussed  in  Section  3.1.1,  the  hydrogen-fed  SOFC-GT  hybrid 
system  can  only  output  70%  of  the  design  power  with  normal 
part-load  strategy.  In  the  present  paper,  we  take  the  system  fed 
by  hydrogen  to  explain  the  three  cases  considering  its  lowest  net 
power  output.  Fig.  12  shows  the  minimum  and  maximum  power 
and  efficiency  curves  of  the  discussed  measure.  Fig.  12  shows  the 
power  output  of  the  hybrid  system  operating  with  different  cases. 
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Fig.  9.  Electrochemical  performance  and  operation  parameters  of  SOFC-GT  hybrid  system  fed  by  hydrogen. 


Case  1:  As  the  turbine  inlet  temperature  is  below  its  design 
value,  auxiliary  firing  upstream  the  turbine  could  significantly 
increase  the  GT  power  output.  An  attempt  has  been  made  by 
adding  auxiliary  fuel  flow  such  that  the  relationship  between 
turbine  inlet  temperature  and  air  flow  equals  the  methane  fuel 
case.  At  the  design  flow  rate,  to  lift  the  TIT  to  1173K,  appro¬ 
priately  15%  of  the  fuel  was  directly  fed  to  the  combustor. 
However,  burning  fuel  generally  increases  heat  production;  thus, 
the  SOFC  fuel  utilization  must  be  decreased  to  keep  the  temper¬ 
ature  constant.  The  simulation  shows  that  with  this  method,  the 
power  can  only  be  increased  insignificantly  to  72.6%  as  shown  in 
Fig.  13. 

Case  2:  To  decrease  heat  recuperation,  part  of  the  exhaust  can  be 
bypassed  around  the  HE1  by  using  the  bypass  valve.  The  simu¬ 
lation  shows  that  this  method  can  increase  the  power  output  at 
the  design  SOFC  operating  temperature  to  a  value  of  84%  of  the 
methane  design  net  power  if  70%  of  the  exhaust  gas  is  bypassed 
as  shown  in  Fig.  13.  However,  net  efficiency  decreased  to  48%, 
caused  by  the  higher  exhaust  temperature.  The  turbine  inlet  tem¬ 
perature  decreased  further,  which  led  to  less  power  contribution 
of  the  gas  turbine. 

Case  3:  If  70%  of  the  exhaust  gas  is  bypassed  and  the  TIT  is 
increased  to  its  original  design  value  of  1 1 73  K  with  auxiliary  fir¬ 


ing,  a  net  power  output  of  appropriately  94%  of  methane  design 
value  can  be  achieved,  as  shown  in  Fig.  13.  Net  efficiency,  how¬ 
ever,  decreased  to  45%. 

Case  3  proved  it  could  increase  the  power  output  almost 
to  its  original,  methane-based  value,  at  the  cost  of  significant 
power  efficiency.  A  promising  control  strategy  for  a  system  mainly 
operating  at  part-load  would  be  to  operate  normally  with  a 
closed  bypass  valve,  and  without  auxiliary  fuel  up  to  a  power 
of  70%,  and  for  a  higher  power  to  gradually  open  the  bypass 
and  introduce  fuel  to  the  combustor  directly.  To  implement  this 
strategy,  the  exhaust  gas  bypass  and  the  auxiliary  fuel  flow 
must  be  integrated  into  the  operating  control  strategy.  Assum¬ 
ing  a  digital  control  unit,  this  does  not  imply  any  hardware 
modification. 

Fig.  14  shows  the  compressor  operation  line  of  different  cases. 
Higher  turbine  inlet  temperature  with  Case  1  and  Case  3  are 
acquired  with  auxiliary  fuel,  which  caused  a  higher  pressure  ratio 
at  a  constant  air  flow  rate.  Therefore,  the  surge  margin  is  decreased 
but  does  not  cross  the  surge  line;  thus,  the  gas  turbine  operation 
is  safe.  Net  power  produced  by  the  gas  turbine  is  higher  than  that 
of  Case  2  and  normal  operation.  In  Case  2,  TIT  is  decreased  rapidly 
due  to  the  bypass  valve,  which  makes  the  compressor  surge  mar- 
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Fig.  10.  Electrochemical  performance  and  operation  parameters  of  SOFC-GT  hybrid  system  fed  by  ethanol. 


Fig.  11.  Gas  turbine  operational  line  with  different  fuels  at  part  load  operations. 


Fig.  12.  Enhancing  power  output  for  the  hydrogen  fueled  system  with  different 
methods. 
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Fig.  13.  Gas  turbine  operation  line  with  different  enhancing  power  cases. 


Fig.  14.  Power  output  of  hybrid  system  operating  with  different  cases. 


gin  increase  at  part-load  operation.  However,  gas  turbine  output 
power  is  thus  reduced  to  a  low  level  as  shown  in  Fig.  13. 

4.  Conclusions 

The  performance  of  an  SOFC-GT  hybrid  system  designed  for 
methane  fuel  operating  with  other  suitable  fuels  was  analyzed. 
New  design  points  were  defined  for  each  fuel.  At  new  design  points, 
the  net  electric  efficiencies  are  62.29,  56.59,  and  60.5,  for  methane, 
hydrogen,  and  ethanol,  respectively. 

The  maximum  output  power  of  the  system  fueled  by  hydrogen 
will  decrease  to  105  kW  (70%  of  the  methane  value)  due  to  lack  of 
cooling  by  stream  reforming.  Due  to  the  absence  of  endothermic 
reforming  reactions,  the  maximum  radial  temperature  gradient  is 
much  higher  than  that  of  the  methane  system. 


When  the  system  is  fueled  by  ethanol,  the  net  output  power  will 
decrease  to  88%  of  the  methane  value  due  to  the  lower  cooling  effect 
of  stream  reforming.  However,  with  increased  heat  recuperation 
from  the  gas  turbine  exhaust  gas,  the  net  efficiency  can  remain  at 
60.5%  at  a  high  level.  This  result  is  based  on  the  assumption  that 
ethanol  vaporization  is  accomplished  by  exhaust  gas  heat.  If  an 
electric  heater  was  used  for  this  purpose,  the  net  power  would  be 
reduced  to  80.02%  of  the  methane  value  and  the  net  efficiency  will 
decrease  to  55.4%. 

To  increase  the  power  output  of  the  hybrid  system  operating 
with  non-designed  fuel,  three  possible  measures  are  introduced  in 
this  paper.  Simulation  results  show  that  for  a  hydrogen-fueled  sys¬ 
tem,  simply  using  the  auxiliary  fuel  to  improve  the  TIT  to  its  design 
value  cannot  increase  the  system  output  power  significantly,  while 
bypassing  gas  exhaust  around  HE1  can  increase  the  net  output 
power  to  84%  of  the  methane  designed  power.  When  both  measures 
are  adopted  at  the  same  time,  the  net  output  power  will  increase  to 
94%  of  methane-designed  power.  However,  the  net  efficiency  will 
decrease  to  45%. 
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